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Introduction 

In  recent  years  renewal  theory  has  become  one  of  the  mo3t  power¬ 
ful  tools  of  the  applied  probabiliato  In  particular,  it  plays  a 
prominent  role  in  the  analysis  of  the  behavior  of  type  I  and  type  II 
particle  counters 0 

There  is  an  Integral  formula  (due  to  R.  Pyke)  for  the  distribution  of 
the  time  between  successive  registrations  with  a  Type  I  counter  under  the 
assumptions  that  the  particles  arrive  according  to  a  general  recurrent 
process  and  that  the  counter  has  a  random  dead  time.  Unfortunately  there 
is  no  such  simple  formula  in  the  case  of  the  Type  II  counter.  In  chapter 
I,  I  give  a  resume  of  the  renewal  theoretic  approach  to  particle  counting 
problems  and  a  shorter  proof  of  Pyke's  formula. 

The  physical  literature  has  been  mainly  concerned  with  the  case 
of  Poisson-arrivals.  The  need  for  considering  general  recurrent  input  is 
genuine,  however,  because  scaling  circuits  and  other  devices  used  by 
physicists  destroy  the  Poissonian  nature  of  the  input  (arrival)  process. 
For  instance,  if  we  count  every  r-th  arrival  in  a  Poisson  process  we 
actually  have  sn  Erlang-r  input  process  (i.e„,  the  time  between  successive 
arrivals  has  an  Erlang-r  distribution). 

In  his  investigations  of  telephone  traffic  fluid  Type  II  particle 
counters  Professor  L.  Tcdc^cs  discussed  both  the  distribution  of  V(t)  = 


the  number  of  registered  particles  at  time  t  and  of  ^(t)  *  the  number 
of  impulses  present  in  the  machine  at  time  t  .  He  formulated  a  recurrence 
system  of  Integral  equations  for  the  binomial  moments  of  the  distribution 
of  ^(t)  and  solved  them  in  the  cace  where  the  dead  time  (impulse  time) 
produced  by  each  particle  is  exponentially  distributed  and  the  input  is 
recurrent,,  Although  the  random  variable  V'(t)  is  of  primary  concern  in 
particle  counting,  ^t)  1®  essential  random  variable  in  the  theory 

of  infinitely  many  server  queueing  systems  as  -£(t)  is  precisely  the  size 
of  the  queue  at  time  t« 

At  the  suggestion  of  Professor  L„  Takac3,  I  attempted  to  utilize  his 
methods  to  solve  the  more  realistic  problem  where  the  dead  time  is  allowed 
to  havo  Erlang-r  distribution  r  >  1,  rather  than  the  exponential.  Unfor¬ 
tunately,  the  equations  become  quite  involved  and  I  succeeded  in  determining 
only  the  Laplace  transforms  of  the  number  of  impulses  present  in  the  counter 
at  time  t  for  Erlang-2  dead  times.  Luckily,  I  discovered  another  class 
of  distributions  which  are  more  tractable  with  respect  to  these  problems, 
even  though  they  are  not  easily  handled.  The  distributions  I  use  are 
wmax~m"  distributions,  the  maximum  of  m  exponential  distributions  each 
with  parameter  p  . 

After  determining  the  number  of  impulses  present  in  the  type  II  counter 
(the  queue  size  problem)  I  turn  to  the  main  problem  of  particle  counting, 
the  behavior  of  v>(t).  Here  is  where  the  "max-m"  distributions  make  the 
problem  manageable.  In  order  to  obtain  the  mean  time  between  consecutive 
registrations  consider  tho  following  bulk-queueing  problem:  Suppose  that 
the  particles  arrive  in  batches  of  size  m,  each  particle  independently 
produces  an  impulse  whose  length  Is  exponentially  distributed  with  parameter 
\l  .  Then  If  we  allow  the  2  stochastic  processes: 


1)  Recurrent  input;  max-m  impulse  time  and 

2)  Recurrent  input;  batch  arrivals  of  size  m  with  each  particle  producing 
an  exponentially  distributed  impulse  to  occur  simultaneously  we  see  that  the 
counter  is  free  at  the  same  time  in  both  processes.  Thus  if  we  let  8(t)  » 
the  number  of  Impulses  present  at  time  t  in  the  second  process,  in 

but  at  any  time  t  when  'jj(t)  =0  8(t)  also  will 

equal  0  and  conversely. 

In  chapter  1  determine  the  binomial  moments  of  the  ergodic  distri¬ 
bution  of  the  imbedded  Markov  Chain  for  the  batch  arrival  problem.  Using 
the  limiting  distribution  of  this  chain  and  Wald’s  Fundamental  Identity  of 
Sequential  Analysis  we  can  determine  the  mean  time  m  between  successive 
registrations  with  max-m  dead  time.  By  the  elementary  renewal  theorem  the 
number  of  "counts"  at  time  t  Is  asymptotically  equal  to  0 

Chapter  5  Is  concerned  with  two  attempts  to  obtain  formulas  for  the 
variance  of  the  time  between  consecutive  registrations.  I  was  unable  to 
obtain  a  formula  in  the  case  of  general  inter-arrival  times,  but  I  give  an 
approximate  result  for  Erlang-r  inter-arrival  distributions  if  r  Is  large. 

My  final  chapter  is  devoted  to  a  modified  version  of  a  problem 
discussed  by  Professor  W.M.  Hirsch  concerning  the  application  of  queueing 
theory  to  missile  defense  systems.  I  treat  the  problem  of  attacking  a 
well  fortified  base  (one  with  infinitely  many  missile  batteries).  Because 
of  the  simplicity  of  the  model  assumed  we  were  able  to  introduce  and  calcu¬ 
late  explicitly  a  loss  function  which  yields  a  reasonable  criteria  for 
evaluating  the  effectiveness  of  certain  attack  and  defense  strategies.  It 
is  hoped  that  this  type  of  loss  function  will  be  "calculable"  In  more 
complex  situations. 


general  6(t)  t  ^(t) 
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Chapter  I  Basic  Renewal  theory  and  the  Type  I  Counter 

The  well  known  Poisson  process  Is  a  model  for  an  Integer  valued 
random  process  {5(t)  j  t  >  0$  which  counts  the  number  of  random  events 
occurring  In  the  time  Interval  (0,t3  .  Usually  the  events  are  represented 
by  the  times  ...  of  their  occurrence.  The  random  variables 

~  *^*1  f  ■*  ^  |  ...  ,  ®  ^  * 0  *  are  called  the 

successive  inter-arrival  times.  More  generally  If  the  time  differences 
are  assumed  to  be  identically  distributed  positive  random  variables  with 
a  common  distribution  function  F(x)  ,  and  if  we  denote  by  the 

number  of  events  which  occur  in  the  time  interval  (0,t3  ,  then  we  say 
that  {©n,  5^1  forms  a  recurrent  (or  renewal)  process. 

Notation;  1)  P£en  —  ** 

oo 

2)  E(e)  =  -C  a  f  xdF(x) 

J0 

3)  V(0)  =<r2  a  f°  (x-x)2  dF(x) 

k)  0(s)  =  e“sxdF (x) 

0 

Definition;  The  renewal  function  (or  mean  value  function)  of  the  renewal 
process  is  defined  by 

m(t)  =  E[$(t)]  =  y  nPfo  ”  n3  • 

n=0 

Thus  m(t)  is  the  expected  number  of  events  occurring  in  the  interval 

(0,t]  o 

The  role  of  the  renewal  function,  both  in  theory  2nd  practice,  can 
hardly  be  overstreased.  According  to  W.L.  Smith  (Ref.  1  page  246^ 


"•••  In  most  applications  of  renewal  theory  a  knowledge  of  the  renewal 
function  m(t)  or  even  a  knowledge  of  lt9  asymptotic  behavior  for  large 
values  of  t,  answers  most  of  the  questions  we  are  likely  to  ask." 

In  fact  knowing  m(t)  we  can  determine  F(x)  from 

5)  £  =  !%• 

The  moat  basic  results  of  renewal  theory  will  now  be  stated,, 

Theorem  1  (Elementary  Renewal  Theorem) 

6)  11m  SliSi  =  i  where  *  =  E(©)  <  oo 

t-O-oo  Z 

and  the  limit  is  interpreted  as  0  if  «<  «  oo  « 

Theorem  2:  If  d  (t)  denotes  the  variance  of  ^  and  if  <T  =  V(0)  <  oo 
then 

7)  11,  HfM.  J£? 

t-^  oo  z 

Theorem  3  (Asymptotic  Normality) 

If  tf-2  <  oo  then  for  all  real  x 

g  4 

8)  11m  pf  - —  <  *}  =  — ±_  fX  a  a“ 

vnn^  "  JZtT- 

Theorem  (Blackwell's  theorem) 

If  the  inter-arrival  time  ©  is  not  a  lattice  random  variable  and 
=  E(8)  <  oo  then  for  any  h  >  0 

9)  11m  Sitekl  ■-  a  i 

t-J>  00  “  * 


Theorem  $  (W.L.  Smith)  If  g(u)  has  bounded  variation  in  the  Interval 
jo,  oo)  and  P(x)  Is  not  a  lattice  distribution  and  Its  mean  •<  <  oo, 
then  we  have 

t  i  co 

11m  \  g(t-u)dm(u)  ■  2  f  g(  )du 
t->ooJ0  J0 

We  now  describe  the  renewal  theoretic  approach  to  the  problems  of 
particle  countingo  Assume  that  particles  arrive  at  a  counter  at  times 
^  a  0  <t1  <T2  <1?  <  ...  <  Tn  <  . . .  where  the  ®n  =  Tn  -  are 

independent,  identically  distributed  positive  random  variables.  Since 
most  counters  have  a  positive  "resolving  time"  not  all  the  particles  that 
arrive  are  counted.  Let  the  subsequence  of  ^  |  denoting  the  arrival 
times  of  the  particles  actually  registered  be  0  =  ^  4  <'£g,<, ,« 

f  y  | 

Again  the  ®n  =  Tn  -  the  times  between  successive  registrations  are 
identically  distributed  positive  random  variables.  Therefore  the  primary 
renewal  process  {©n , generates  a  secondary  renewal  process  l®n*  s> 
where  v*t  *  the  number  of  particles  counted  In  time  (0,t}„  We  denote 
by  R(x)  the  probability  P£®n'  <  xj  . 

In  order  to  ascertain  the  asymptotic  behavior  of  the  number  of  recorded 
particles  p(t) ,  it  suffices  to  determine  the  common  mean  m  and  variance 

p 

c T  of  the  distribution  function  R(x),  Once  wo  know  these  values,  then, 

by  the  basic  renewal  theorems  given  earlier 3 

ita  .  1  lta  =  jd 

t->  00  z  m  t— >00  t 

The  mechanism  of  the  counter  used  determines  how  the  subsequence 

of  "recorded  events"  is  selected  from  the  primary  sequence  of  events 


We  say  that  the  counter  is  free  at  time  fc,  if  it  is  in  condition 
to  register  a  particle  arriving  at  that  time ,  otherwise  the  counter  i3 
said  to  be  locked.  The  two  standard  counter  mechanism  are  the  type  I 
and  type  II  counters.  In  a  type  I  counter.  If  a  particle  arrives  when 
the  counter  is  free  then  the  counter  Is  locked  for  a  random  time  X  called 
the  dead  time  or  holding  time  or  Impulse  time  of  the  counter.  Particles 
arriving  during  the  time  %  are  not  counted  and  have  no  effect  on  the 
counter’s  operation.  The  type  II  counter  differs  from  the  typo  I  counter 
in  that  every  arriving  particle  locks  tho  counter  for  a  random  locking 
time  %.  .  Therefore  a  particle  now  is  registered  if  and  only  if  at  the 
time  of  it3  arrival  ati  the  dead  times  produced  by  the  previous  arrivals 
have  expired.  With  a  type  I  counter  a  particle  is  registered  if  at  the 
time  of  Its  arrivals  all  the  dead  time  produced  by  the  previous  registered 
particles  have  expired.  It  is  assumed,  of  course,  that  tho  dead  times 
produced  by  the  arriving  particles  (typell  counter)  or  by  the  registered 
particles  (type  I  counter)  are  identically  distributed  independent  random 
variables. 

When  the  particles  arrive  according  to  a  Poisson  process,  i.e„# 

F(x)  =  P [@n  <  xj  -  1  -  e  '  the  following  results  are  known,  (see  Takacs 
[5*J  and  Smith  [4]  ) . 

For  a  type  I  counter  m  =  -  <j-^  _  jLJlA— XiZi. 

*  A2 


where  X  Is  the  dead  time  random  variable  P[Z<  x]  =  H(x) 
Thus 

E[V(t)l  .  A  V[p(tfl  \  1+A2V(*) 

lim  11  — *•  -  rrgarwrt  J-im  — “P — »  =  A  ftr-ViAf’ 


oo 


rrmr 


For  a  type  II  counter: 


i  „*E  M 


while 


6AeKI 

“T 


4  "pUi 


<?P 


1  -  H(y)  dy  -  lj  dt  + 


2,XlsfC}.2  Bfia 

— Ji — 


In  order  to  present  our  shorter  proof  of  Pyke's  and  Ma3®.qni$t  »s 
formulas  for  the  typo  I  counter  we  must  give  one  more  definition* 
Definition:  If  wo  havo  the  first  registration  at  Q  ■  0,  an  event  of  the 
pi  umary  sequence  is  registerable  If  it  arrives  after  the  dead  time  Xq 
produced  by  the  first  particle „ 

Theorem  6  (Pyke 1  s  Formula):  If  Pjl^  *^h-i  —  *  F(x)  and  the  ®n  =£T^ 

are  mutually  Independent  and  If  the  dead  times  produced  by  the  n-th 
particle  (if  It  is  registered)  are  distributed  as  H(x)  and  the  £&n$  ave 
mutually  independent  and  also  are  independent  of  the  sequence  iv  • 

Then  0  *  =  has  the  distribution 

z 

p[©n  —  ZJ  *  R(z)  =  P-  *  P(z-u)3  dm(u) 

where 

oo 

s  pn(x)s 
n*L 

Proof:  Clearly  we  can  assume  that  there  is  a  registration  at  time  0 

I0e.  T0  =  #t0'  =  0  .  Then 

R(z)  *  P  [there  is  a  registration  In  (0,z)3 

=  P (there  is  a  regiaterable  event  in  (0,z3) 

=  P[the  last  event  in  (0,zl  is  registerablej 

CO  «  3>I  /  H(u-)£i  -  p(z-u)}  dPV(n). 
n*T  0  n 

(\JU)  mf  H(u-)  [l  -  P(z-u)]  dM(u)  o 

Jo 


The  reasoning  is  as  follows:  the  probability  that  the  n-th  event  is 
the  last  one  occuring  in  (0,z3  6111(1  18  regiaterable  is  by  the  total  pro- 

i  * 

bability  theorem  H(u-)[l  -  F(z-u)j  dPn(u)  (because  the  n-th  event 
occurs  at  time  u,  0<usas).  As  the  last  event  in  (0,j£|  can  be  the  n-th 
for  n  *  0,1,2,  . ..  again  by  the  total  probability  theorem  we  obtain  that 


op  z 

R(z)  *  Z>  f  H(u~)ll  *•  F(z-u)}dF  (u) 
n^T  J0  n 

*  H(u-)  Ql  -  P(z-u)3dM(u) 

Theorem  7»  (Malmqulst’a  Formula).  Under  the  same  general  assumptions  of 
Pyke’s  theorem,  except  that  now  we  assume  that  the  counter  is  locked  for 
a  constant  time  d  when  it  registers  a  particle,  we  have 

(3)  R(z)  *jO  s  <d  ^ 

l?(z)  -  F(d)  +f  {f(z-x)  -  P(d-x)}dM(x)  . 

J0 

Proof : 

R(z)  *  P|Jbhere  is  a  regi3terable  event  in  (0,jQ) 

**  P Qihere  is  a  first  regiaterable  event  in  (0,z)J 
Here  the  first  regiaterable  event  is  th»  first  one  to  arrive  after  time 
d  »  The  first  arrival  is  the  first  regiaterable  if  it  arrives  in  (d,z3« 
The  probability  of  this  is  just  F(z)  -  F(d)„  If  n  >  2,  the  n-th  particle 
is  the  first  regiaterable  if  the  (n-l)-st  particle  arrives  at  time  x  , 
x  <  d  and  the  n-th  arrives  between  d  and  z  By  the  total  probability 
theorem  the  probability  that  the  n-th  (n  >  2)  is  first  regiaterable  is 


therefore 


d 


Applying  the  total  probability  theorem  again  (on  n)  we  obtain 

oo  £ 

R(z)  *  F(z)  -  F(d)  +  21  f  fP(s-x)  -  F(d-x)}dP  ,  (x) 

nS?  Jo  n”1 

«  F(z)  -  F(d)  +  j{P(t-x)  -  F(d-x)}  dM(x)  . 

0 
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Chapter  II 


O/Eg/oo 


Introduction: 

Consider  a  type  II  counter.  The  particles  arrive  at  the  instants 
Tg,  ...  where  n  ■  1,2,...  are  identically  distributed 

positive  random  variables  with  distribution  function  F(x).  Let 

1)  <P(s)  =  f  e"sxdF(x) 

°0 


We  shall  assume  that  F(x)  is  not  a  lattice  distribution  and  that 
«<  <  oo  ,  Recall  that  in  a  type  II  counter  every  arriving  particle 
produces  an  impulse  but  only  those  particles  arriving  whan  the  counter  is 
free  will  be  registered.  Let  'Xn  denote  the  impulse  time  or  dead  time 
produced  by  the  n-th  arriving  particle.  We  suppose  that  the  n  «  1*2,. 

are  identically  distributed  mutually  independent  positive  random  variables 
which  are  also  Independent  of  «  In  the  present  chapter  we  suppose 

that  the  ^  have  an  Erlang-2  distribution. 


1+)  H(x)  =  P  <  x)  *  1  -  -pxe* 


■px 


and  we  shall  derive  a  formula  for  the  Laplace  transform  of  the  r-th 
binomial  moment  of  the  number  of  impulses  present  In  the  counter  at  time  t 
This  random  variable  will  be  called  ^p(t)  .  Although  theoretically  we 
can  invert  the  binomial  moments  and  obtain  the  exact  distribution  of^t) 
the  fomulas  obtained  are  too  complicated  for  practical  use  in  this 


manner. 


The  random  variable  ^p(t)  is  not  very  interesting  from  the  point 
of  view  of  particle  counting  but  is  the  object  of  main  interest  in  the 
theory  of  infinitely  many  server  queueing  systems  since  ^(t)  is  now 
the  queue  size  of  the  system  G/Eg/  oo  • 

Notation:  P^(t )  =  PC^(t)  =  k3 

The  r-th  binomial  moment  cf  «^(t)  is  given  by 


oo 


Br(t)  »  21 0  Pv(t) 


By  Jordan's  inversion  formula 

00 


Pv(t)  =  T  M)r‘k  {I)  Bjt) 


00 


B. 


(s)  =  f  e"8^  B  (t)  Re(s)  >  0 

Wo  r 


r 


is  the  Laplace  transform  of  tho  r-th  binomial  moment  of  ^(t)„ 


oo 


OP 


=  f  xdH(x)  =  f  fl  -  H(x)Jdx  will  always  be 

*»0  J0 

assumed  to  exist  and  be  finite  for  all  dead  time  distributions  used  in 
this  paper. 


§  2 .  A  Review  of  the  Results  of  Takucs 

In  this  section  we  shall  assume  that  0) 
Consider  the  generating  function 


0  and  Pjh^  <  x]  *  F{x) 


(2,1)  G(t,a)  = 


oo 


Pk(t)  zl 


Theorem  1:  The  generating  function  G(t,z)  satisfies  the  following 


Integral  equation 


2.2.  G(t,z)  ■  £  -  F(t)]  +  j£  G(t-x,z){z  +  (1-z)  H(t-xi)}  dF(x) 

Proof? 

By  the  theorem  of  total  probability 

P0(t)  *  1  -  F(t)  +J^  H(t-x)  PQ(t-x)  dF(x) 

since  we  have  no  impulses  present  in  the  machine  if  either  there  is  no 
arrival  up  to  time  t  (the  probability  of  this  event  is  1  -  F(t))  or 
there  is  a  first  particle  arriving  at  x,  0  <  x  <  t  ;  the  impulse  it 
produces  expires  by  time  t  (the  probability  of  this  is  Just  H(t-x))  and 
the  process  "renewed”  at  x  has  no  impulses  present  at  time  t  0?o(t-x)3  . 
Similarly  we  obtain  for  k  *  1,2,  ... 

(2o4)  Pk(t)  *JoiPk(t-x)  H(t-x)  +  Pk-1(t-x)  (l  -  H(t-x)J$dF(x) 

Multiplying  the  equations  (2.I4.)  by  z^  and  adding  over  k  =  0,1,2,  ... 
we  obtain  (2.2)  as  desired. 

Theorem  2.  The  binomial  moments  Br(t)  exist  for  all  t  and  csn  be  determined 
from  the  following  recurrence  formulas: 

Bq ( t )  =1  and 
.t 

(2.5)  Br(t)  *  J  Br-1(t-x)  jj  *•  H(t-xQ  dm(x)  r  =  1,2,  ... 

00 

where  m(x)  *  Fn(x)  and  Fn(x)  denotes  the  n-th  -  fold  convolution  of 

F(x)  with  itself.  Further 


Proof:  Clearly  Br(t) 


1 


r  =  0  ^1  ^2  f  * » • 


Since  B0(t)  *  1  upon  differentiating  equation  (2,2)  r  times  with 
respect  to  z  and  evaluating  at  z  *  1  we  obtain 


jj 

Bit)  =  f  B_(t-x)  dF(x)  +  j*  B  , (t-x) jl  -  H(t-x)]dF(x) 
r  wq  v  v0 


This  a  linear  integral  equation  of  the  Volterra  type  for  the  unknown 
B  (t ) o  The  solution  is  well  known  (it  ta  obtainable  by  taking  Laplace 

P 

transforms)  to  be 


V4’  =  1 


'0 


(l  -  H(t-x)3  dm(x)  . 


To  prove  (206)  we  must  show  that  the  generating  function  Q(t,z)  is 

Cr 

analytic  at  z  »  1,  It  suffices  to  show  that  Bp(t)  <  ^  for  some 
constant  C.  By  (2,5)  we  can  write 

Br^  c  <  t  ^  "  H(tg-t^)3  *  * 4  [1  “  H(t-'tr)J  dm(t^)o» 

Let  b  be  a  fixed  positive  number  and  let  k(x)  *  H(x-h)  ,  Since 


m(t+h)  -  m(t)  <  1  +  m(h)  for  all  t  >  0  we 


easily  obtain 


V«  J-L  -  k(*lO  -  &  - 

+  *  •  •+xr5t+h 


<  00 
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<  p^gM|r  where  J>*  ^  xdH(x) 


thus  setting  C  *  .(h+»  we  have  proved  that 

Cr(t,z)  is  analytic  at  z  =  1  and  thus  it  is  permissible  to  invert 
oo  . 

Br(t)  *  (£)  and  we  obtain 

-  g  (-Dr-k  (£)  Br(t)  . 

Before  proceeding  to  discuss  the  limiting  behavior  we  review  some 
results  from  the  theory  of  functions  of  bounded  variation,, 

Lemma  l:  The  product  of  two  functions  each  of  which  ic  of  bounded  variation 
is  also  of  bounded  variation. 

Lemma  2:  The  integral  of  a  function  of  bounded  variation  on  a  finite  interval 
[Vb]  *3  also  of  bounded  variation  on  that  interval 0 

Theorem  3 •  If  T*  xdH(x)  <  oo  xdP(x)  and  F(x)  is  not  a 

-  J0  J0 

lattice  distribution  then  the  limiting  distribution 


lim  Vi'ijit)  *=  k}  ~  Pj*1  k  =  0,1,2,  exists 

t->  co  1  * 

and  is  independent  of  the  initial  distribution  and  we  have 

«•«  \*  “  %  <-Dr"k  <S>  Br* 

where  Bp*  is  the  r-th  binomial  moment  of  {P^*}  and  can  be  determined  in 
the  following  way: 


(2.8)  Bq*  »  1  and  Bp*  ■  i  J"  Br-1(t)  |l  -  H(t)]  dt 
Proof: 

First  we  shall  prove  the  theorem  in  the  particular  case  that  f0)  " 0  • 

As  BQ(t)  -  1  clearly  B *  1  and  by  (2.5)  we  have 

t  - 

B1(t)  *  j  [l  -  H(t-xQ  da(x) 

As  1  -  H(x)  is  monotone  non-increasing  vie  may  apply  Smith's  theorem 
to  deduce  that  exists  and  equals 

v  *  j  r  &  -  H<t)]  dt  *=  . 

By  lemma  1  B^(t)  is  of  bounded  variation  on  every  finite  interval  and 

as  the  limit  exists  B^(t)  is  bounded  on  the  entire  line,  we  can 

show  (2.8)  by  induction.  If  we  assume  that 

lim  B_  .  (t)  =  B  exists,  then  by  Smith's  theorem 
t->oo  r“1  r“1 

,t 

applied  to  Br(t)  =  j  B^-^t-x)  (l  -  H(t-x)3  dm(x) 
we  can  conclude  that 

lim  B  (t)  a  B  *  exists  and  is  given  by  (2.8)  0 
t— >  00 

I* 

Since  Br*(t)  <  —y  for  t  >  0 

(2.9)  lim  G(t,z)  *  0*(t) 
t-»>oo 


exists  and  we  have 


The  aeries  (2,10)  is  convergent  for  all  z.  By  (2,8)  we  see  that 


(2,11)  3*(z)  •  1  -  P°  o(t,z)  (l  -  H(tQ  dt, 

*  Jo  "" 

Now  G*(l)  =  1  and  according  to  the  continuity  theorem  for  generating 
functions  it  follows  that  the  limiting  probabilities 


lim  P£»(t)  =  k]  =  Pk*  k  =  0,1,2,, 
t— >oo  1 


exist  and  that 


G*(z)  * 


00 


*  _k 


P,„  z 


Finally  from  (2,10) 


«  *  _  1  dV*(z) 

\  *  *r  ~zr^ 


•  fl  (-l)1'11  (£)  B 

s=0  r3£  k  r 


and  thus  the  Bp*  (r  =  0 ,1 ,2 ,  , , , ) 

are  indeed,  the  binomial  moments  of  the  distribution  iPk*J  ,  This 

completes  the  proof  for  the  initial  state  0,  If  we  consider  an  arbitrary 

initial  state  then  the  only  difference  is  that  ^(t)  is  replaced  by 

'ijft-'t^)  +  €(t)  where  is  the  random  arrival  time  of  the  first  new 

customer,  ^(t)  has  the  same  distribution  as  our  ^(t)  and 

lim  P 0 (t )  s  0 ]  »  1  ,  (€(t)  is  the  number  of  the  original  customer 
t— >oo 

oo 

remaining  in  the  queue  at  time  t.  Clearly  l  xdH(x)  <  oo  implies  that 

0 

€(t)  ->0  as  t  — >oo)  .  Consequently,  ij(t-T^)  +  €(t)  has  the  same 
asymptotic  distribution  as  the  special  case  »  0  that  we  considered. 

Thus  the  proof  of  the  main  theorem  is  now  complete. 
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j  3<-  The  explicit  Solution  In  an  Operational  Form 

The  fundamental  recurrence  between  B  (t)  and  B„  .  (t )  la 

Br(t)  a  B^tt-x)  £l  -  H(t-xO  . 

When  H(x)  is  Erlang-2  this  becomes 
(3ol)  Br(t>  a  f  Br-1(t-x)  fo”^***5  4-  p,(t-x)  ©‘'^(t-xij 

Upon  taking  Laplaco  transoms  we  obtain  the  recurrence  equation  for  the 
corresponding  transforms: 

(3.2)  Br(s)  =  CBr.i(»+»)  -  Br^  (s+n)]  . 

As  the  units  in  which  we  measure  tine  are  completely  arbitrary  we 
may  assume  that  p  =  1  without  loss  of  generality.,  In  order  to  make  the 
discussion  as  simple  as  possible  vie  introduce  the  following  notations  and 
operations o 


E:  shifts  r  ~>r  +  1  e0g»  B^(a)  =  EP0(a)  ;  py(s)  =  E^Bgts) 

FS  shifts  the  argument  3  — >  s  +  1  e„g<,  Fd(s)  =  d(s-KL)  ;  F^dCs)  =  d(s+k) 

Ds  is  the  differentiation  operator . 

It  is  Important  to  note  that  while  "0,E  and  F  commute  with  each  other  , 
they  do  not  commute  with  functions  of  3,  Our  equation  (3»2)  now  becomes 
the  operational  equation: 

(3.3s 


Er  pQ(s)  =  d(l  -DjFE1’”1  Pq(s)  given  P0(s)  =  >s 


Theorem  3-1 


(3c3)  E1’  «  JZ  «(*i  .... ,  a JdD^FdTi*2?2 d„,.Dar-1li,r*1d(l‘D)arPr 


(*1  *•  •  °  , 


ap>l 


r 


i  “i  -  3  *»  "  3  -  ”  '  1 

where 

(3olj.)  C  f  o  •  e  “  ( '*1 


Before  giving  the  proof  we  gf.ve  some  preliminary  explanation.,  To 


find  E 


such  thai. 


,r .  1 


•/e  )  ss  p  (a)  we  E.ust  first  enumerr.te  all  the  r-tuples  (a-,  .  ,,8^.) 

a,  -  v 


.1  _ 

t  1)  a  ,  >  1.  2)  ®;l  <  J  foi>  3  <  *’  -  1  3)  ^ 


For  each  r-~tuple  consider 


(3.5)  cO^joou  ,ap)  .d(s)  D^dCu-KL)  D^dU^)  ...  D  ^d  .'a+(r-l) )  Cl-D)  ’rd(s*r) 


where  c(a.  ,0»,,aJ  is  a  constant  depending  on  the  arrangement  c  ho  son  end 
(1-D  )R  i.j  to  be  expanded  ir.  the  binomial  formula.  B ( 3  is  equal  bc>  the 


uum  of  expressions  of  the  form  of  (3»5>)  . 

It  Is  into re a ting  to  identify  the  combinatorial  meaning  of  the  terms 
in  the  product  of  t be  coefficient  c  (a^  »» .  *  ,2.p)  „  For  choosing  a^  we 
sire  given  one  ball  jj  wo  can  take  it  (taka  1 )  or  leave  it  (taka  Oi  For 
choosing  we  are  given  a  new  ball;  If  *  0  we  choose  ag  balls 

from  two  balls  i.e.  ag  =  0,1,  or  ?.  while  if  *  1  ag  **  0  or  1. 
Similarly,  we  choose  balls  from  j  balls  minus  the  number  of  balls 


The  product 


croaon  p..-avJous»ly  ..  .l.o » »  choose  a,’  from  j  • 

C  4=2,  \ '  1 

V  ""i  (  4  •».  i  J 

/  equals  the  lumber  of  ways  wo  can  put  r*  .India  fcJn^uiahal.lo 


i 


.1 


J-JL  V  * 

balls  .Into  j?  colls  l:a  such  a  way  that  tho  number  of  balls  in  the  first 
j  colls  is  always  <  j  » 

Proof  rf  ths  thoorom; 

T!- ...  mc'hccl  proof  i.s’  by  induction.  A**  tho  formula  holds  for 


■?  . i • r  *  epci « i~d ) ■»  <roFd < i«d ) p*1 


.sfiimi-  tint  It  :1s  titio  for  B1  ***  and  show  it  holds  for  .IS**  „ 


-•  d(l«P)P  '^Z 


a 


3*Fd., .  “d  (l-.D)2Fr‘  J 


^w,aJdD  - 


( a;i  •  •  •• "  #*V-2 


.f  “•  ,r: 
NT' 


u,  *  a  « 
T  * 


a.  <  -j 

to  x  - 


a  >  l 


for  till  i  «:  v 


O 


;?o  x;Ua‘fc  :1m  happening  look  rvi  the  :!n  F; 

a. 


ni 

c  (a^  . ,  ...  ,p.)  dC  Fd  * „  D 


:’^"2a(w>)8  p5’"1 


oott  b\':. tv  i ,  o 


C  (q.^  j*  >  *:•  -  * 


'v 


a ‘ 


Best  Available  Copy 


0(0, 

*1.  ••• 

fX  <►)  J  0  s>  (I  *£l  „  />  O 

f.  .‘■‘-'V  ' 

c(l. 

:11! 

:a^2,  a) 

0(0, 

V 

1  ;<  fi.p  *  S>  3  «  0-^.,  ^  1  1 

o  (C, 

V 

B<*-  X  ^  *)  !t  U  '  U  »  A 

c.  ;>  y>%*tL 

o(0. 

al  •■’ 

j>d  (i  A  ^  i<^  ji  p * 

•??r. 


(we  omitted  the  operations  since  the  a's  in  the  e-term  determine  the 
corresponding  operation).  By  inverting  the  above  system  and  remembering 


that 

*  1  op  0  we  see 

that 

(3.6) 

c(l#  •  •  •  •  #  i 

a)  =  -  c(a2>**»# 

while 

(3.7) 

C  V  0  9  &2  y  o  *  o  f 

,  a)  «  +  c (aj *  ar«l * 

®  c (ag-l t  S3 }  « *  o  f  a) 


*  c  ( fip }  }  aj }  a ) 

o 

-  c(&2#  •••  i  *  at*«*l*^ , 

Notice  that  if  «  2  c(ag»  a3-l,...)  will  be  zero  as  (g)  *  0  and 
In  general  (^)  =  0  If  k  >  n,  With  this  convention  it  is  easy  to  see  that 
(3»6)  holds  l.So 


(»  1) 


21  a 


c  (1  $  &/)  i«  •  a  $  &)  ®  •  c  (&2  jm>  » o  i  y  n*“X )  sine© 

r-1 


1  OCl2) ...  c 


r-2 

r-2  -  >  a 


a2  x  a3 


‘r-1 


‘).  a 


V3 


7^ same  product] 


To  see  that  (3.7)  holds  for  our  formula  we  use  the  elementary  relation 


Thus  : 
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(- 1) 


2 


a2  -  fk  -  a3  -  «2 

a3  /  V  *4 


Proceeding  in  this  way  (i.e.,  decomps® 


term  next  etc.) 


we  obtain 


r-1 


(al 


90  o  o  |C 


r-l  t 


a)  =  (-  1) 
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Hence 

c(0#  a2,...,  a^,  a)  *  -  c(a2  -  1»  a^,  ...»  a) 

—  c  (a2  ,  a^  —  X  |«  nf  a) 


—  c  (a2  •  a^ i  • « t  f  4^ ^  —  1|  a) 

—  c  ( &2 *  a3 j • • • i  5  a  —  X ) 


as  desired  and  we’ve  shown  that  the  c's  fit  the  relations  (3.6)  and  (3.7) 
thus  if  (3.3)  holds  for  E1*"1,  it  holds  for  E1*. ,  Therefore,  as  (3.3)  is 
correct  for  r  «  1,2,  by  the  induction  hypothesis  E?  is  given  by  (3.3) 
for  all  r  . 


4.  The  Limiting  Distribution  of  ^(t)  » 

By  standard  Tauberian  theory  (see  e.g,  Doetsch  PeQ )  if  lim  B  (t) 

Woo  r 

exists  then  lim  B  (t)  *  lim  sp  (s)  •  Since  we  know  (f 2)  that 


Woo  *  s->5  ”‘'r 


B  «  lim  B  (t)  exists  and  that  {b  *?  for  r  =  0,1,2, 
t->oo  r 

binomial  moments  of  the  limiting  distribution  P.  *  lim 

*  Woo 


...  are  the 

Pfy(t)  =  k)3 


k  =  0,1,2,  ...  we  can  apply  the  Tauberian  theorem  referred  to.  We  now 
compute  Bj*  ,  B2*  and  B^*  . 

Lemma  1:  lim  a€ts)  =  ~  where  «<  *  \  xdP(x)  <  00 

s  -*0  T-?sT  *  JO 


1  -  *8  +  0(8^) 
*<8  +  O(s^) 


Proof: 
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11m 

s-*0 


r 


The  first  three  binomial  moments  of  the  limiting  distribution  are 


1) 


Proof: 


11m  Sp,(s)  *  11m  T.ffTaT 

8-^0  1  3-^0  A  ^3; 


s+2  _  2 

7 


2) 


Proof:  B2(s)  *  dFd(l-D)2P2(1/fl)  -  dDPdd-DjP2  X/a 


B2(a)  =  I%sT  °  l'-^TT  C1  -  20  +  °23  -  r-^fsT1 


s  iUl  .  « -jSB&tlLjL-  nJL  +  — — — +  2  j .  jfliL  i—£—)  rd  ^  "I 

i“fls)  l-fTs+IT  ^a+2  (s+2)^  (s+2)^  l-fia)'!-^'  '•s+2  TZIpY^  j 


(s+2)‘ 


therefore 


Similarly  It  ean  be  shown  that 


3>  B3«.icr m  -i m  -  %  rM-  »-3 

+  rfiJj)  '  m  (t£%)'  <2)3 
+  ^C  Ui^)^1)^)  (2)+T^T7  WiS^)  (2)3 

In  the  special  case  of  Poisson  arrivals  £(s)  **  Va+s  and  we  obtain 

B^*-  2X  B2*=  bX2  B,*. 

which  agree  with  earlier  results  of  Tallies  who  showed  that 

**  *s  lim  P^ft)  *  e**^  ^r-)*-  where  J°  =  mean  impulse  time  *  2 

W  no 

if  the  locking  time  has  an  Erlang-2  distribution,,  Since  the  binomial 

r 

moments  of  the  Poisson  distribution  of  density  V  are  g.Iven  by  Bp  *=  v  /rl 
our  formulas  coincide  with  those  of  Takacs, 
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Chapter  III 


Q/nax-o/oo 


1  Introduction 

In  this  chapter  we  again  determine  the  Laplace  transforms  of  the 
binomial  moments  of  ^ (t )  *  the  number  of  Impulses  present  In  the  counter 
at  time  t  *  the  queue  size  of  the  queueing  system  O/max-m/oo  .  Thus 
we  suppose  that  particles  arrive  at  a  type  II  counter  according  to  a 
recurrent  process  where  the  inter-arrival  times  -  ^n-1^  51  *  1»2,»»«> 

are  independent  positive,  non-lattice  random  variaDles. 

Let  F(x)  =  P[rn  -rn-1  <x] 

oo 

-c  =  f  xdF(x)  <  00 

<(s)  *  T  e~sxdF(x) 

J0 

Finally  we  assume  that  the  locking  time  ^  produced  by  the  n-th  partiole 
is  the  maximum  of  m  independent  identically  distributed  exponential  random 
variables  l<>e<>, 

H(x)  =  P[*n  <  d  -  &  -  •'>“  3"  ■*  ^  (-  1)J  (  5  )  .-1’*3 

It  will  be  seen  that  the  recurrence  equations  for  the  Laplace  transforms 
of  the  binomial  moments  will  be  easier  to  solve  for  the  general  max-m 
distribution  than  for  the  case  of  Erlang-2  dead  time  distribution  dealt 
with  in  the  previous  chapter,, 


2 •  The  Transient  Behavior  of  "jj(t) 

In  1958  CO  Professor  L.  Takacs  derived  the  following  recurrence 
equations  for  the  r-th  binomial  moment  at  time  t. 
t 

(2*1)  Br(t)  BjJt  -  x)  £l  -  H(t-x)J  dm(x) 


where  H(x)  is  the  deadwtime  distribution.  In  our  situation  the  equation 
to  be  solved  becomes 


(2.2) 


m  }  e-n(t-x)jj 


Letting  0  (s)  =  f  e 
r  «/0 

we  obtain  from  (2*2) 


Bp(t)  dt  be  the  Laplace  transform  of  Bp(t) 
the  equation 


(2.3)  0rU)  -  (-  l)3*1  ( 

for  the  Laplace  transform  of  the  r-th  binomial  moment*  It  Is  to  be 
noticed  that  no  derivative  appears  In  (2*3)  In  contrast  to  equation  (3*2) 
of  Chapter  II  and  thus  we  have  a  pure  m-tk  order  difference  equation  to 
solve*  As  BQ(t)  *  1  for  all  t  PQ(s)  =  for  all  s.  For  convenience 

we  shall  set  n  =  1,  Y^fs}  3  1111(31  \  3  d(a+k)  •  Then  (2.3) 

becomes 

Pr(.)  -  4(.)  (-  l)3*1  (  “  )  Pp.1(.+J)} 


(2.U) 
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Theorem:  he  solution  of  (2.4)  is 


(2.5)  Pr(s) 


(-1)^ 

StK 


d°  <3tJl  dui+u2*  du1+...+up_1  i«i  ^  uj 


over  all  partitions 
(u^,  ... ,Uy)  of  k 
such  that 

;>  "  Uj  =  k  J  1  <  Uj  <  m 


Proof: 


If  r  =  0  Pn(s)  »  Vs.  If  r  =  1  from  (2.4) 


(2.6) 


=  40  §  <-l>  ,’1(  mi  )  S+T 


Now  there  is  one  and  only  one  partition  of  J  satisfying  the  conditions 
of  (2.5),  namely,  j  itself 0  Therefore  the  coefficient  of  *3 

<-  I)**1  dQ(  ^) )  as  desired  and  thus  formula  (2.5)  gives  the  solution 

if  r  «  1  .  We  proceed  by  induction.  Assume  that  the  theorem  is  true 
for  all  moments  up  to  r  -  1  .  In  particular 


m(r-l)  ,  1 jk+r-1 

for  «ir 


r-l  m 

doW+M2  •••  dW«°‘+ur-2  i£  <ul 

over  all  partitions 

(^1  ,u2  9  •  •  •  *  ^-l  ^  of  k 
r-l 

such  that  !j>  u^  *  k 
and  1  <  u  j  <  a 


Since  0p(s)  *  d(s)  it  (s+J it  la  fairly  clear  that 

Pp(s)  is  given  by  our  formula  because  the  shift  a  +  k  — a  +  k  +  J 

has  many  partitions ,  each  of  whioh  comes  from  a  chosen  from 

1,  and  another  partition  Ug»  •••  ,  Up  which  is  a  partition  of 

k  +  J  -  u-.  o  The  proper  multiplier  in  each  case  is  (  ®  )  or  "fr  (  ®  ) 

1  ui  i=2  ui 

Thus.,  all  we  have  done  is  to  introduce  new  partitions  of  the  numbers  r  to 

mr  the  first  of  which  is  j  « 

In  detail: 


(2.7) 

m  «  i  «  m(r-l) ,  1  \k+r-l 


dodj\.j  •  Vu^-  •  •  V  j  K 

over  all  partitions 
<«!*...  .Up^)  of  k 

such  that 
r-1 

Uji  *  k  1  <  u^  <  m 


m  m(r-l)  /  i)k+r-j  f  __ 

s+k+J  V  ^  j*‘*du1+u2+. 

\  over  all  partitions 

(v^^Ug,... ,Up-;L)  of  k 

such  that 
r-1 


_  r-1 

(  m  )  ~TT  l  m 
*  +ur-2+J '  J  1=1  * 


-  JJ- 


Lot  ¥  =  k  +  J  t  then  ¥  varies  from  r  to  mr  and  our  sum  becomes 


s+y 


d„d..  d..  ...  d.. .  jt  <  5 ) 


"•  V— ■ %.i 


1=1 


over  all  partitions  of  y 
such  that 

^  “l 


as  desired.  In  words,  the  j  of  the  first  sot  (2.7)  becomes  the  first 
element  Uj^  of  the  partition  of  y  into  r  parts. 


/  3o  The  Ergodlc  Behavior  of  ^(t)  . 

Although  the  formulas  are  rather  cumbersome  we  can  determine  the 
first  few  moments  of  the  limiting  distribution  and  can  use  the  standard 
Tauberian  theory  to  obtain  a  general  formula  Br*  the  r-th  binomial 
moment  of  the  limiting  distribution  11a  P[  ^(t)  =  kf)  (see 

chapter  II  /U.  Since  Br*  exists  it  is  given  by  a*  =  lim  sp  (s)  . 

r  3  0  r 

Hence 


(3ol)  Br* 


1 

7 


mr 


(-i)k+P 


>  8  8 

«1  ul+u2 

all  partitions 

..Up)  of  k 


ul+u2+> *  »+ur-i 


such  that 


k 


1  <  Uj  <  m 


i 


where  5^  * 


To  illustrate  the  use  of  (3.1)  we  compute  Bg  for  the  case 
whore  H(x)  is  max  2. 


B.. 


=,  1  A-  i^l.Lkl2.  ^3 


k 


5  C 

U-,  U-.  -'U, 
X  ..  (. 


(  ®  )  (  m  ) 

«»  tu 


over  all  partitions 

^ul,u2^  of  k 
•such  that 

u^  +  ru,  =  k  1  <  <  2 


the  partitions  (u^jUg)  of  k 

!\  «  n  n  j, 

n  n  n  n 


k  * 


2  aro  «  1 ,  tu  ®  I 


3  n 

4  " 


«  1#  Up  *  2  u^  «  2,  u^ 


=  1 


*1 


-  2,  u.,. 


therefore  3, 


*  1 


{- 


2  i 


(f)  (i> 


{2.  ,2> 

l2; 


<!><!> 

“T~ 


=  *  ^  °1  ”  ®2 


If 

0,. 


^(o)  a  Va+s  »  that  is  if  the  arrivals  form  a  Poisson  Process  then 
■  Vk  ana  we  obtain 


v 

B„  *  /V 


r^-4A3  * 


(4  A)' 

t r 


as  desired 


Of  course  **  where  S  is  mean  nlocking  time'*  and  «<  is  mean  service 

time  for  general  distributions. 

References:  Same  as  Chapter  II0 


Chapter  TjT  :  A  related  Batch-Arrival  Queueing  problem  and  the  deter¬ 
mination  of  the  mean  time  between  consecutive  registrations, 

£  1  Introduction;  So  far  we  have  concentrated  our  attention  on  the 
random  variable  -^(t)  s  the  number  of  impulses  present  in  the  machine 
at  time  t  .  Although  the  queue  size  ^(t),  is  of  fundamental  import¬ 
ance  in  the  theory  of  infinitely  many  server  queues  it  is  not  usually 
of  primary  interest  in  the  theory  of  particle  counters,  'As  we  remarked 
in  the  Introduction,  the  random  variable  central  to  the  theory  of 
Counters  is  the  "time  between  successive  registrations'1  which  we  denoted 
by  R(x)  o 

In  order  to  ascertain  the  mean  of  R(x)  when  the  dead  time  has 
a  max-m  distribution  we  introduce  a  nev;  Stochastic  Process  which  may 
also  be  of  interest  In  the  theory  of  bulk  queues.  We  shall  assume  that 
particles  arrive  In  batches  of  size  m  at  the  type  II  counter  each 
particle  producing  an  exponentially  distributed  impulse  Independently  of 
the  other  particles.  Let  8(t)  be  the  number  of  impulses  present  In 
the  counter  in  this  new  process.  In  general,  8(t)  is  not  equal  to 
(t) ,  but  B(t)  and  4£(t)  will  be  zero  simultaneously.  Thus  the 
time  between  two  arrivals  finding  the  counter  free  (i.e0.  R(x))  will  be 
the  same  In  both  cases.  We  shall  see  that  it  is  relatively  easy  to 
find  the  mean  of  R(x)  in  the  8(t)  process  because  we  can  utilize  a 
l-dimensional  imbedded  Markov  Chain, 

f  2.  The  General  theory  for  the  Bulk-Arrival  Model 


In  this  section  we  shall  derive  the  recurrence  equations  for  the 
batch  arrival  system  where  the  particles  arrive  at  the  type  II  counter 


(or  customers  at  an  infinitely  many  server  queueing  system)  in  batohes 
of  size  m  according  to  a  recurrent  process •  The  inter-arrival 
distribution  F(x)  Is  assumed  to  be  non-lattice  and  positive.  The 
"dead  time"  produced  by  any  particle  (or  service  time  of  a  customer) 
is  assumed  to  have  the  positive  distribution  function  H(x)  .  Let 
J  be  the  sequence  of  arrival  points  of  the  groups  of  particles  and 
let  %  be  the  dead  time  produced  by  a  single  particle.  We  assume  that 
P[Tn  <  x]  =  F(x)  for  all  n  and  that  the  t'f  are 

mutually  independent.  Also  the  individual  dead  times  produced  are 
mutually  independent  and  also  independent  of  the  iT*n^  system.  Let  £(t) 
denote  the  number  of  impulses  present  in  the  counter  at  time  t  (equivalently 
£(t)  denotes  the  queue  size  at  time  t).  For  simplicity  we  assume 
that  S(0)  =  0  and  P[*f^  <  x][  =  F(x).  If  we  let  Pj£(t)  =  P[fl(t)  =  kj  we 
can  derive  an  integral  equation  satisfied  by 


G(t  ,z) 


Pk(t)zk  . 


Theorem  1:  The  generating  function  G(t,z)  satisfies  the  following 
integral  equation 

_  Jfc 

(2.1)  G(t,z)  =  Jl-F(tO  +  f  G(t-x,z)[z  +  (1-z )  H(t-x)3“dF(x) 

Proof:  By  the  theorem  of  total  probability 

(2.2)  PQ(t)  ■  (i  -  F(t)]  *f0  mpQ(t-x)  dF(x) 


(2.3)  Pk(t) 


“  )lH(t-xQm"J  [1  -  H(t-xflJ  PkfcJ(t-x)dF(x) 


1  lie  1*3  fore 
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oo 


(2,14.)  G(t,z)  »2^Pk(t)zk 


Il-P(t)|  +  J^  ^  (“)»k[l.H(t.x)]\l(t.x)f  J.Pk-j(t-x)ak*JdP(x) 


Interchanging  summation  yields 


G(t,z)  =  [l~F(t)]  +  j*  ^3?  o(“)zJ  jl-H(t-x)]J(H(t-x)],n”Jpk-J(t-x)zk’JdP(x) 


or 

^  BX 

G(t,z)  =  fl  -  P(t)]  +  f  G(t-x,z)  fz  +  (1-z)  H(t-x)l  dP(x) 

J0 


(4*60! 


If  we  differentiate  G(t,z)  and  evaluate  at  z  =  1  then  by 
exactly  the  same  method  as  used  In  chapter  2  in  the  section  summarizing 
Takaca *  results  we  can  prove: 

Theorem  2:  The  binomial  moments  exist  and  can  be  obtained  by  the 
recurrence  formula. 


(2.5)  Br(t)  =  ®  ^  Br_j  (t-x)  [l-H(t-x)J^  dP(x) 


Cor.  The  mean  of  the  ergodic  distribution  Is  simply  where 

op 

T  *  mean  dead  time  =  f  xdH(x)  <  00  and  •<  *  mean  inter-arrival 

J0 

00 

time  *  f  xdP(x)  <00  . 

Jo 


Proof:  As  BQ(t)  *  1, 


Bl(t)  *J*  {^(t-x)  +  m[l-H(t-xO]  dP(x) 
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thus  taking  Laplace  transforms  yields 

B^s)  *  Br(s)  f(s)  +  m  f(s)  L  [1-bQ 

where  L[l- h]  is  the  Laplace  transform  of  1  -  H(x)  * 

Hence  P^(s)  =  m  Y'-ffs)  **  15-  “  *0  and  ^  inversion 

(l  -  H(t-xf}  dm(x) 

and  F  (x)  denotes  the  n-fold  convolution 
n 

of  F(x)  o 

Applying  Smith’s  renewal  theorem  yields 
(2,6)  B1*  *  .  q.e.d. 

The  major  problem  in  using  theorem  2  for  practical  purposes  is 
the  presence  of  the  J  in  the  exponent  of  (l  -  H(t-x)3  in  formula 
(2o!>).  This  difficulty  will  be  high-lighted  when  we  specialize  to  the 
case  of  exponential  dead  times  in  the  next  section.. 

ho  The  Special  Case  of  Exponential  Dead  Time  and  the  Mean  registration 
time  for  the  type  II  counter  with  Max-m  dead  time* 

In  this  section  H(x)  *  1  -  e***  so  that  (2<,E>)  becomes 
(3.1)  Br(t)  «  7  )  Br_J  (t-x)  (t"X,dP(x)  . 


h.™  -  Jo 
00 

where  m(x)  =  Fn(x) 


Upon  taking  Laplace  transforms  we  obtain  the  recurrence  equation 
for  the  Laplace  transforms  of  the  r-th  binomial  moments.  Namely, 

(3»2)  pr(s)  =  T^sT  ^  {  J  }  Pr-J  (a+^) 

It  13  to  be  noted  that  this  equation  will  in  general  be  an  m-th  order 
difference  equation  with  variable  coefficients  and  will  be  difficult 
to  solve  unless  m  **  1  , 

We  now  proceed  to  discuss  the  Imbedded  Markov  Chain  in  the  case 
of  exponential  dead-time  0  We  shall  again  see  that  the  technical 
difficulties  in  this  approach  are  due  to  the  fact  that  the  m  will 
appear  as  an  exponent  in  the  integral  equation  determining  the  generating 
function  of  the  queue  size  at  arrivals. 

If  we  let  6n  *  6  (Tn  -  0)  be  the  number  of  Impulses  present  Just 
before  the  n-th  batch  arrives  the  sequence  \J>n\  forms  a  Markov  Chain 
with  transition  probabilities 

(3.3)  P(6ntr  =  nten  »  3)  *  Pjk  ■=  f  (3J“>  a'1"11  (i-.n«)3«>-*  . 

If  we  can  find  one  solution  to  the  system  of  equations 

(3.U)  Pk"  ^  ?3p3k 

such  that  {Pk\  Is  a  probability  distribution  then  by  Poster's  theorem 
QQ  the  chain  will  be  ergodlc  and  1 P^\  will  be  Its  actual  unique 
limiting  distribution,  independent  of  Initial  conditions.  Proceeding 


00 


as  usual  we  introduce  the  generating  functions  U(s)  *  «  Since 


*  ■  SL  fo  <J~> 


J-ha-k 


dP(x) 


Multiplying  by  zK  and  adding  over  k  we  obtain 

oo  go  J+m  J+m-k 

U(z)  =2”  l  >  P,  (JJ“)  e"^  zk  (l-e^)  dP(x) 

jSff  J0  J  K 


or 


(3o5)  TT(z)  =  (1  -  ^  +  se-^)m  TJ (1  -  e~^  +  ze^)  dp(x) 


oo 


Recalling  that  B  =  > _  v  /  x. 

i*  1 


(*)  Pv  is  also  equal  to 


1  drtT(z) 

F*  tor  »=i 


and  differentiating  with  respect  to  z  r  times  in 


the  given  equation  we  conclude  that  the  binomial  moments  of  the  limiting 
distribution  satisfy  the  difference  equation 


(3.6) 


Bp  ■  tw  ^  (P!j)  Bj  =  nvM)  ^  >  Vj 


If  we  can  solve  (3.6)  for  the  binomial  moments  then  we  can  find  the 

p 

probabilities  by  Jordanfs  inversion  formula.  In  particular  ®/«< 

is  the  mean  recurrence  time  of  the  time  between  batches  arriving  and 
finding  the  counter  free  (or  all  servers  free  in  the  case  of  the 


r0 

Infinitely  many  server  queue).  However,  this  /•<  Is  also  the  mean  time 
between  successive  registrations  In  the  ^(t)  process  as  i^(t)  *  0  Iff, 
6 (t )  *  0  ,  This  Is,  of  course,  the  reason  we  considered  the  batch 
arrival  problem. 

Notice  that  if  we  again  set  \l  =  1  d^  * 


(3,6)  becomes 


(3.7) 


(3  *  d 

V  V 


{^r  (",B^  lf  r-n 


(3.7) 


B  B  d 
r  r 


m 


<  1  > 


lf  r  >  m 


and 


1 1 

mTT*[ 


B„  = 


€2 


Theorem!:  The  general  solution  to  (3.7)  is  given  by 


(3.8)  B„ 


\ 


d  * 

**•  S+u2+*“+"k  ill 


(  “  ) 


u. 


all  ordered  partitions  of  r 

into  r,r-l,...,l  parts. 

i.e.  for  each  k  <  r  consider 
~  k 

partitions  (u^. .. 

1  <  Uj  <  m 


»  r 


J , 


Proof:  By  definition  BQ  *  1  «  B^  *  mdj^  which  is  also  given  by  our 
formula  since  there  is  only  one  partition  of  1  namely  u^  *  1  • 

For  Bg  there  are  exactly  two  partitions  of  two  (1,1)  and  (2)  so 
our  formula  yields 

(  ®  J^d^d^  +  ^2  ^d2  agreeing  with  a  "brute 

force"  computation*  The  general  proof  is  by  induction*  Suppose  the 
theorem  is  true  for  all  r  <  r  -  1  we  must  show  it  is  true  for  B„  . 

P 


If  r 


-■  Br  =  "  >Vj* 


therefore 


B. 


r 


k 

El  ddu  ...  d  1 T(  “  > 

«i  ui  *,,  uk  i=i  *1 

of  all  ordered  partitions 

of  r-J  into  k  parts  k  =  1,2,... »r-j 
k 

such  that  ^Uj,  =  r  «  j 
1  <  Uj  <  m 


k 

— rp/  m 

i=l  V^^ij+Ug  •••  du1+. 


■+uk 


all  ordered  partitions 

of  r  into  k  parts  k  =  1,,..  ,r 
k 

'*1 


^’r 


where  the  empty  bracket  la  to  be  read  as  one.  Notloe  that  the  (  ®  ) 

term  is  Just  the  last  part  of  the  partition  of  r  Into  k  parts  l.e.  , 
Ujj  *  J  and  thus  we  see  that  the  proof  If  r  <  m  Is  complete.  The 
proof  for  r  >  m  Is  similar.  As  the  formula  (3.8)  Is  true  for 
r  *  m  by  the  previous  proof  we  again  use  Induction. 


If  r  >  m 


21  II  (  n.  )du-du-+u- 

1=1  *1  "l  V*2 

all  ordered  partitions 

of  r-j  into  k  part3 

k  =  1,...,  r-j  such  that 
k 


*+uk 


since  the  (  ®  )  term  again  becomes  tho  last  part  J  of  the  total  partition 
of  r  and  combines  with  all  partitions  of  r-j  to  give  a  partition  of 
r  into  k  +  1  parts.  Note  that  for  some  k,  especially  small  ones, 
there  may  be  no  relevant  partitions  as  the  partition  of  r  may  have 
Uj1 8  >  m  and  (™)  =0  if  u^  >  m  .  Actually  it  is  not  surprising 

that  (3»7)  and(3.7’)  have  the  same  solution  as  they  are  identical  if  we 
make  the  convention  that  13^  =  0  If  k>u. 

We  now  are  in  position  to  give  the  formula  for  the  mean  of  R(x). 


Theorem  U.:  The  mean  time  between  successive  registrations  in  the  type  II 
Particle  Counter  where  the  arrivals  form  a  recurrent  process  i.e. 


<  x)  »  P(x)  -t  ■  xdF(x)  <  oo 


are  max-m  distributed  is 

A 

V 


0 

where 


and  the  impulse  times 


P 


0  * 


(-1)% 


where  the  By  ere  given  by  (3*9)  . 

Proof:  By  the  fundamental  theorem  of  x*ecurrent  events  applied  to 
Markov  Chains  the  mean  number  of  steps  In  the  chain  before  a  return 
to  Eq  starting  from  Eq  Is  ^/Pq  *  However  the  length  of  time  between 
successive  steps  of  the  chain  is  distributed  as  F(x)  .  If  we  regard 
the  R(x)  distribution  to  be  the  sum  of  N  F(x)»s  where  N  *  the 
number  of  steps  in  the  chain  (or  number  particles  that  arrive)  until 
an  arriving  particle  finds  the  counter  free  we  can  apply  Wald's 
Fundamental  Identity  to  conclude  that 

oo 

E(R)  »  — as  *<  =  f  xdF(x)  <  oo  . 

*0  Jo 

Notice  that  the  event  In  «s  n)  depends  only  on  the  first  n«l  F(x) 
variables  and  thus  the  use  of  Wald's  identity  Is  Justified,  q.e.d 

00 

(3,9)  Since  PQ  =  ^  (-1)% 

rr 

and  as  B  <  — .  for  some  constant  C  the  series  for  P q  converges. 


The  usefulness  of  our  approach  depends  on  the  rapidity  of  the  convergence, 
not  just  because  we  want  to  add  up  as  few  terms  as  possible  but  also 
because  the  Br’s  are  difficult  and  tedious  to  compute  for  large  r* 
Nevertheless,  by  the  introduction  of  an  auxilliary  stochastic  process, 
we  have  been  able  tc  give  a  precise  formula  "V^q  for  the  mean  time 
between  successive  registrations . 

One  standard  device  used  in  particle  counting  is  a  scaling 
circuit  which  lets  only  every  r-th  particle  through  to  the  counter,, 

This,  mathematically  speaking,  transforms  the  usual  Poisson  input  into 
Erlang-r  input.  In  this  case 


d(s)  - 


(Va+s )r 

l-(AA+a)r 


and 


B0»l 
*1  = 


m(^A+l)r 

i-(AA+Dr 


B 


2 


m2  rk+ur. 

i-(AA+Dr 


LxA*z_j-  + 

i-(AA+2)r 


- 

l-(A/A+2^r 


If  A  «=  (i  =  1  m  =  2  r~6!j.  then  Bg  <  . 

The  usefulness  of  (3.9)  therefore  depends  on  the  relative  sizes 
of  m  and  r  .  In  order  to  approximate  a  true  dead  time  distribution 
a  large  m  will  be  needed  as 


46- 


the  mean  of  a  max-m  distribution  is 

m  i 

the  variance"  "  "  "  x/ j2 

and  most  dead  time  distributions  will  be  "almost  constant"  and  thus 
the  mean  must  be  much  larger  than  the  variance. 

§  Using  the  ^(t)  Process  to  Investigate  the  5(t)  process  when 
arrivals  are  Poisson 

So  far  we  have  simply  used  the  "batch  arrivals"  process  as  a 
tool  to  determine  the  mean  time  between  consecutive  registrations  for 
the  type  II  counter  with  general  input  and  max-m  Impulse  time.  The 
present  section  Is  devoted  to  a  minor  reversal  of  this  procedure. 

Suppose  we  really  were  interested  in  the  infinitely  many  server 
queueing  process  where  arrivals  obey  a  Poisson  law  with  parameter  A 
service  time  is  exponential  with  parameter  p.  and  the  arrivals  are  in 
batches  of  m  .  The  standard  approach  Qf|  yields  the  following  set 
of  difference  equations  for  the  limiting  distribution  of  the  queue 
size, 

if  k  <  m  -  1  (k  +  DpPj^  =  (\  +■  kp,)Pk 

if  (k+i)pPk+1  =  ( X  +kp.)  -  pk-m 

The  usual  method  is  to  solve  these  equations  in  terms  of  P0 

oo  u 

and  then  normalize  using  the  condition  ^  P^  **  1  ,  In  our  case 


It  is  difficult  to  obtain  the  general  solution  but  PQ  is  known 
to  be  equal  to  e”  where  *  *  mean  of  the  max-m  distribution. 
The  Justification  for  this  assertion  is  that  if  we  consider  the 
two  processes  in  operation  simultaneously  we  are  in  0  in  the  $ 
process  if  and  only  if  we  are  in  state  0  of  the  ^  process 
G/max  m/oo  .  By  TaKacs'  result  [l]  PQ  *  e’^  .  In  fact,  Takacs 
shows  that  the  limiting  queue  size  for  an  infinitely  many  server 
queue  with  Poisson  arrivals  and  general  service  distribution  H(x) 
is  Poisson  with  parameter  9  *  |  xdH(x)  <  #  , 
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Chapter  V*  Attempts  to  Determine  the  Variance  of  the 
"Time  between  Consecutive  Registrations." 

fla  Introduction: 

In  order  to  obtain  the  asymptotic  distribution  of  , 
the  number  of  registered  particles  at  time  t  ,  we  need  not  only 
the  mean  of  R(x)  but  also  the  variance  of  R(x).  Although  we  have 
not  been  able  to  obtain  a  formula  for  the  variance  of  R(x)  ,  we  can 
give  an  approximation  to  the  variance  for  Erlang-r  inter-arrival 
distributions  if  r  is  large. 

Despite  the  fact  that  we  were  unsuccessful  in  our  attempt 
to  discover  the  variance  of  R(x)  «  it  is  informative  to  review  the 
approaches  used  and  to  see  the  difficulties  involved.  The  first 
method  uses  Wald' 3  Fundamental  Identity.  Let  us  consider  the  S  pro¬ 
cess.  When  the  first  batch  of  arrivals  comes  we  keep  on  sampling 
from  the  inter-arrival  distribution  until  an  arriving  batch  finds 
the  counter  free.  If  N  *  the  number  of  arrivals  betwean  the  regist¬ 
rations  including  the  2nd  registered  particle  Q-.e.  N  =  n  means  that 
after  registration  1,  rv»l  particles  wsre  not  counted  but  the  n-th 
was]  ,  then  N  =  the  recurrence  time  of  the  state  EQ  of  the  imbedded 
chair  •  If  ©  stands  for  the  inter-arrival  distribution  and 

Sr  »  ®i  +  Sg  +  •  •  o  +  ©jj  then  Sr  has  the  distribution  R(x).  We 
used  this  idea  previously  in  showing  that 

E(R)  »  *7*0  aa  "■*  a™*  E<N)  “  Vp0  • 


Assuming  that  both  N  and  9  have  finite  second  moments  we  can  differ' 


entiate  Wald’s  Identity  (see  Harris  [l]  3  obtain  the  following 
formula  for  : 

E[S/J  *  2E(NSjj)E(©'  -  E[N2JE2(®)  +  E(H)V(9) 

The  main  difficulty  encountered  -in  this  approach  is  the  determination 

p 

of  E(NSjj)  o  The  quantity  E(N  )  is  also  hard  to  compute;  nevertheless 

p 

we  can  give  a  formula  for  a  (N)  and  thus  for  E(N“)  0  Because 

02(N)  i3  very  complicated  this  method  will  not  be  of  practical 
use  ever*  if  we  could  solve  for  e(ns^~|.  It  is  of  interest  to  realize 
that  the  random  variable  H  can  be  defined  to  be  the  number  of  steps 
between  successive  transitions  EQ  — >Eot  in  the  chain  ° 

In  the  second  part  of  the  chapter  we  outline  an  approach 
using  integral  equations,  also  developed  by  Takacs.  Unfortunately, 
this  second  approach  yields  less  than  the  first,, 

§ 2o  An  Approach  via  Wald’s  Fundamental  Identity  and  Recurrent 

Events  o 

A3  already  indicated,  the  present  approach  Is  an  attempt  to 
compute  the  second  moment  of  R(x)  by  use  of  the  formula 

(2ol)  E[R2}  *  E[Sn23  *  2E[NSN]E(e)  -  e[n2]E2(®)  +  E(N)  V(0) 

We  shall  first  determine  O^CN)  »  E(N2)  »  E2(N)  ,  the  variance  of 
the  number  of  steps  between  consecutive  transitions  EQ  in  the 


Imbedded  chain  . 

Recall  that  pjk  *  (^)e”likx  (1  -  •~lJX)^‘m"kdP(x)  . 

Let  fp]/n^}  be  the  distribution  of  the  number  of  impulses  present 
in  the  8  process  Just  before  the  n-th  arrival.  Starting  from 
the  initial  distribution  {Pk^i  the  distributions  {Pk^}  can 
be  determined  recursively  by  the  formulas „ 

(2d)  Pv(n+1)  *  l?”  p  P.(n)  n  =1,2,3,  ... 

K  3SE-m  rJk  J 

Consider  the  binomial  moments  of  the  distribution  of  the  queue  size 
at  the  n-th  step 

(2.2)  Br(n)  =  E  [  (pn)$  *  g  (k)  Pk(n)  . 

Using  the  elementary  result  that  the  binomial  moments  of  the 
Binomial  (Bernoulli)  distribution  Qk  *  (k)plc(l-p)n“k  are  just 

Br  *  we  can  Prove; 

Theorem:  BQV =1  n  =  1,2,  ...  and 

(2.3)  Br(n+1)  =  <?p  (  ®  )Br_k*n)  where  B^n)  =0  ifV>  0 


Proof:  If  we  let  ®n  ■  Tn  *  Tn_^ 


9 


as  usual ,  then  conditional  upon 


en  *  x  and  »n  =  J 


n 


,3+m)#-na 


because  under  the  given  conditions  5^^  is  just  a  Binomial  variable 
with  parameters  j  +  m  and  o”^3c  thus 


B  (n+1) 
r 


£ 


(3«,Pj(n) 


e 


AS  (J«) 


331 


(  ^  )  (m) 
'r«k'  'k; 


Br(n+1)=  <Pr  ±(?){±  <rWn^ 


B 


(n+1)  _  <p  »  (n) 

1  "  ^r  «&*V  r-k 


Q  e  e  o  d  a 


it  (n )  « 

If  n  ~ >*oo  and  B  »  lim  B  v  '  then  the  B  satisfy 

r  n~>oo 

equations  previously  derived  for  the  ergodic  distribution  (see  Chapter 
TV  equation  (3®6)).  Starting  from  B  ^  r  »  1,2,  ....  the  binomial 
moments  of  the  number  of  impulses  present  at  the  first  arrival  we  can 
use  (2o3)  to  obtain  the  binomial  moments  of  the  number  of  impulses 
present  just  before  the  n-th  arrival  (remember  we  are  working  with  the 
5(t)  process)  »  If  0(0)  *  i  and  ^  =  x  ,  0^  has  a  Bernoulli 


distribution  with  parameters  i  and  e"^  and  thus: 


(Zok)  Bj/Hejf*3)  -  (J)^,  r  *  0,1,2 


>  •  •  • 


We  proceed  to  determine  the  generating  function  for  the  r-th 
binomial  moment  of  {6^$,  Lot  Bp(w)  *  By ^ w^  0 

Lemma:  Suppose  6(0)  »  i  then  the  Bp(w)  satisfy  the  difference 

equation: 


<2-5>  V">  ■  ~r!^  { 


m 


*k*Br-k*w)} 


Proof: 


Bp  ^n+1  ^  *  Cj. 


-®-  ,m>R  (n) 
VBr-k 


Multiplying  by  and  summing  over  n  we  obtain 


oo 


B.(n+1)«n*1 

X* 


00 


m 


5:  <?)Erin,wn 


or 


Br(«>  - 1 ;  > v  - 


a 


hence 


btM 


m 


S>wfl 


q.e  «dt 


For  the  practical  purpose  of  particle  counting  theory  we 


-53- 


we  nay  assume  that  initially  we  ape  in  state  Eq  i.e.  6(0)  ■  0  , 
Thus  the  equations  to  be  solved  are 


(2o6) 


Br(w)  -  Vv)l 


B0(K>  -  iSw 

where  dr(w)  =  and  d0(w)  *  . 

‘r 

Theorem:  The  solutions  of  the  equations  (2«6)  for  the  determination 

of  the  generating  functions  of  the  transient  behavior  of  the  binomial 
moments  of  {,6^  are  given  by: 


(2o7) 

Bp(w) 


<*  )dn(w)d„ 


u 


all  ordered  partitions 

of  r  into  k  parts ,  k  «  1 , . , . ,  r 
k 

such  that  Uj  *  r  and  1  <  u^  <  m 


The  proof  is  omitted  as  it  is  exactly  the  same  as  the  one  given  for 
the  imbedded  chain  except  that  d^,  is  now  replaced  by  dp(w)  0  This 
accounts  for  the  presence  of  dQ(w)  *  in  the  formula  (2»7)  « 

As  BQ  s  1,  dQ  *  1  and  so  dQ  failed  to  appear  in  (3o8)  of  chapter  IV  . 
Since  Br(w)  is  analytic  in  a  neighborhood  of  the  origin. 


P0(v)  * 


in  a  neighborhood  of  the  origin*  In  theory,  therefore,  we  can  use 
Cauchy's  Integral  formula  for  derivatives  to  determine  the 

*P0(n)5  n  *  1,2,  , *  Let  us  again  consider  tha*  »w»ry  step  In 

the  chain  is  one  trial  in  a  recurrent  event  scheme  where  the  recurrent 

event  £  Is  "we  are  In  state  0  at  the  n-th  trial"  (or  step  of 

the  chain)  .  The  rs?idom  vai*iable  N  which  equals  the  number  of  trials 

between  successive  occurrences  of  £  and  which  equals  the  number  of 

steps  between  consecutive  transitions  En  — ^  E  in  the  chain  f5  } 

u  in  n 

Is,  of  course,  the  object  of  our  discussion.  We  already  know  that 
E(N)  *  •  By  a  problem  in  Feller's  text  [2"J  we  can  also 

determine  the  variance  of  II  „  Specifically 


(2*8) 


(P, 


(n) 


n 


-  po>  = 


ar^(N)  -  1/PQ  +  (Vp0)2 
- - 


where  Pq  is  the  Pq  of  the  argodic  distribution  and  is  given  in 
section  3  of  Chapter  17  .  Another  expression  for  the  right  side  of 
(2.8)  Is 


(2.9) 


6-2(N)  -  Vp0  +  (Vp0)2 
P^/Pq)2 


Pq)** 


« 


lim 

w->l 


(-1,r  V*’  -  <  ^ 


2 

Unfortunately,  the  determination  of  (N)  In  this  manner 
requires  far  more  computation  than  the  determination  of  Pq  because 

<y2(N)  depends  on  the  sum  (Pq^  •  Pq)  0  Also  even  if  Fq^ 


la  "close  to"  PQ 


many  terms  will  be  needed  to  assure  that 

00  /„) 

is  close  to  (PqV  '  -  Pq).  In  theory,  never 


theless,  we  have  determined  both  the  mean  and  variance  of  the 
recurrence  time  of  the  O-state  In  the  Imbedded  Markov  Chain  of  the 
6(t)  process.  Although  we  were  able  to  deduce  that  the  mean  time 
between  successive  registrations  in  both  the  6(t)  and  the  ^(t) 
processes  is  equal  to  *VPq  ,  we  have  been  unable  to  discover  an 
expression  for  the  variance  of  the  "time  between  successive  regist¬ 
rations."  The  difficulty  arises  from  the  facv  that  Wald's 
Fundamental  Identity  of  Sequential  Analysis  yields  an  intractable 
expression  for  the  second  moment  of  the  sum  of  a  random  number  of 
random  variables  when  the  number  chosen  depends  on  a  sequential 
stopping  rule.  Under  suitable  conditions  (see  Harris  jVJ )  we  have 


(2.10)  E^2]  a  2E(NSn)E(Q)  -  E(N2)E2(6)  +  E(N)V(e) 


where  N  a  the  number  of  random  variables  ©n  ,  all  of  which  are 
independent  and  identically  distributed  as  F(x),  chosen.  The 
term  E(NSK)  is  well  known  to  be  the  "troublemaker." 


In  an  important  special  case  we  can,  however,  approximate 


56- 


Suppose  we  have  a  type  II  counter  with  max-m  Impulse  time  distribution 
and  we  are  counting  particles  arriving  according  to  a  Poisson  Process. 
If  the  particles  pass  through  a  scalar  which  lets  only  every  r-th 
particle  through  then  the  input  to  the  counter  Is  an  Erlang-r  process. 
If  r  Is  large  it  will  almost  be  a  constant  .  now  calculate 

E(NSjj)  if  the  inter-arrival  distribution  is  a  constant. 

Lemma;  If  tho  Input  process  Is  a  constant  i.e.  the  time  between 
consecutive  arrivals  is  .<  >  0.  Then 

(2.11)  sTnSjj!  «  E(N2)-<  . 


Proofs  Under  the  hypotheses  of  the  lemma  PpjN=nJ  m  Pp(H  «  n |Sn^n*Q  thus 


E(NSN)  a  >  n(-<n)Pp(N  «  n,^  =  n»0 
oo 

“  A- 

n*l 


n2-<  Pp(U  *  n)  *  ^E(N2) 


When  the  arrivals  have  an  Erlang-r  Inter-arrival  distribution, 
therefore,  we  suggest  approximating  E(RSN)  by  «<E(N2)  in  equation 
(2.10).  As  E(6)  3  rA  -  and  V(0)  a*  r^2  as  0  i3  now  Erlang-r 
we  obtain  from  (2.10) 

Var(Sjj)  as  eTS^]  -  E2[Sp]  »  2-tE(N2)E(2)  -  E(N2)E2(9) 

+  E(N)V(2)  -  E2(SIJ) 

*  r^2E(N2)  +  E(N2)  +  E(N2)  +  E(N)  ^  -  (*/? Q)Z 

At 


r'ry 


Thus  If  V(0)  ■  r^2  1*  small  (thus  how  large  r  must  be  for 

this  approximation  to  be  used  depends  on  A  )  we  can  say  that 
approximately  d"2(R) ,  the  variance  of  the  time  between  consecutive 
registrations  is: 


(2ol2) 


(Si  * 


^<T2(N) 


Since  the  times  £'t’n  j  of  registrations  of  a  particle  form 
a  recurrent  process,  with  P{Tn^  •  <  xj  *  R(x)  and  if  we 

denote  by  ,  the  number  of  registrations  by  time  t  then 
Feller’s  Central  Limit  theorem  for  recurrent  events  (theorem  3  of 
chapter  I)  asserts  that  if  CT  (R)  <  oo  then  the  distribution  of 
Vt  is  as5rmptotically  normal  i.e0 


(2a13) 


where  p.  a  B(R)  C_^s<T2{R), 

Applying  this  result  to  the  type  II  counter  with  Erlang-r  inter¬ 
arrival  times  and  max-m  dead  time  distributions  the  number  of 

registrations  by  time  t  has  approximately  the  asymptotic  normal 
distribution  with  parameters, 

H  «  n(R)  =*y-  • 

0*2  a®2 (R)  a  «-2(n)  +  -i - ■£- 

X2  Fo  x2 


*  * 


In  the  paragraph  above  we  have  blithely  written  /  Pq  and 
or^(N)  m  It  must  be  remembered  that  even  in  the  case  of  Erlang-r 
input  they  are  not  at  all  easy  to  compute, 

£  2.  The  registration  time  Distribution  in  the  Case  of  Constant 

Inter-arrival  times. 

So  far  w©  have  devoted  our  attention  to  the  case  where 
F(x)  Is  non-lattice.  In  considering  the  variance  of  the  time 
between  successive  registrations  for  Erlang-r  input,  we  used  constant 
arrivals  to  approximate  Erlang-r  Input  and  were  able  to  obtain  an 
approximation  for  O'  (R)  and  thus  give  an  approximation  for  the 
asymptotic  distribution  of  the  number  of  registrations.  When  the 
Inter-arrival  distributions  Is  constant  (•<)  It  is  easy  to  derive  the 
exact  distribution  of  the  time  R  between  successive  registrations. 
VJe  shall  again  work  with  the  8  process ,  the  arriving  particle 
produces  m  Independent  exponentially  distributed  (parameter  p) 
impulses. 

Theorem:  P(R  a  k*<^[  =  7  -(p*)^018 

^1*^2^ *  *  * ^k* 
all  ( *•  •  • » Jjj) 

such  that  <  m  <  2m 

•••  1  < 
and  +. , a  km 


Proofs  In  order  that  the  k-th  arriving  particle  after  a  regist¬ 
ration  be  the  first  one  registered  there  must  be  at  least  one 
impulse  present  at  the  times  when  particles  1,2,.. .  ,  k-1  arrive. 
Therefore,  at  time  •<  no  more  than  (m-1)  of  the  original  impulses 
can  have  ended  i.e.  impulses  can  expire  by  time  •<  where 

*  0,1,..  o,  m-1  but  not  m  .  In  the  next  time  interval  (^s  2«T[ 

^2  impulses  and  but  if  the  second  arrival  i3  not  to  be  registered 
$ 2  <  m~j^  +  m  or  +  Jg  <  2m  »  Similarly  we  see  that  in  the 
i-th  (jkk)  Interval  of  length  *<3  Jj(i-l)*<,I«0j  impulses  and 

satisfies  +  •  •  •  +  Jj_  <  «  Finally  if  the  k-th 

particle  Is  the  first  to  b©  registered,  =  the  number  of  Impulses 
expiring  in  (_(k-l)^,  k*J  is  equal  to  *  km  -  +...+  J^-l^  ° 

Therefore  s 

P[R  =  1*0  =  al 

all  ( •  1 1  j^.)  such  th&t 

<  3^  ***  jg  <  2ru^##oy  < 

3^  +  o «  q  +  3^  ~  km 

2:  I  rr  q.e.a. 

1=1  Ji* 

all  such  that 

<m,..o,  +•••+  <  (k-l)m 

+•  • .+  jjj.  -  Inn 

To  the  best  of  the  author’s  knowledge  there  is  no  known  closed 


form  for  this  sum 
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Chapter  0B: 


The  Integral  Equation  Approach 


Although  we  have  been  able  to  obtain  the  mean  of  the  "time 
between  successive  registrations"  assuming  that  the  impulse  time 
obeyed  a  max-m  law,  we  could  not  obtain  the  variance  since  it 
depends  on  eJns^]  .  In  1957,  Professor  Takacs  suggested  a  new 
approach  which  lead3  to  integral  equations  for  the  mean  and 
variance  of  the  times  between  consecutive  registrations;  conditional 
on  the  impulse  time  produced  by  the  previous  registrant.  After 
summarizing  Takacrf  methods  we  solve  the  integral  equations  for 
a  max-2  distribution.  The  situation,  even  for  max-2  dead  time 
Is  surprisingly  complicated. 

§ 2.  A  Summary  of  Takacs’  Results 

We  assume,  as  usual,  that  particles  arrive  at  a  co’inter 
at  times  ^1*  ^2*  °*°  •  ®ac^  on®  produces  an  impulse  of 

duration  {X^  where  pt^n  ~  =  and  the  [9^  are 

identically  distributed  and  independent  of  one  another  and  of  the 
&n\  sequence.  We  denote  the  subsequence  of  n  =  0,1,2,... 

of  those  particles  actually  registered  by  n  =  0,1,  ...  . 

When  the  times  tr„  -  Vi>  n  *  1,2,  ...  between  successive 
arrivals  are  Identically  distributed  a3  P(x)  then  the  time 
differences  {T*  -  between  successive  registrations  will 

also  be  identically  distributed  random  variables  say  R(x)  .  If 
we  let 


p  *  T  xdP(x) 
J0 


f  (x-p)2dP(x) 

Jo 


then  our  objective  is  to  find 


(2.1)  A  -  m[t£  -  *  Jo  xdR(x) 

(2.2)  B2  =  D2jT^  -  1 (x-A)2dR(x) 

p 

Once  we've  found  A  and  B  we  can  apply  Feller's  Central  Limit 

theorem  of  renewal  theory  to  derive  the  asymptotic  distribution  of 

^  =  the  number  of  registrations  occurring  in  the  time  interval 

2 

(0,tQ,  provided  of  course,  that  A  and  B  are  finite.  To 
2 

determine  A  and  B  we  introduce  the  following  conditional  expect¬ 
ations  o 

(2.3)  A(j)  =  m[T[  |7(0  =  £| 

(2 .k)  B2(y)  =  !5«o  =  rZI 

Knowing  these,  by  the  theorem  of  total  expectation  we  have 

(£.5)  A  *  f°  A(y)dH(y) 

J0 

and 

o  00  «s  oo  ^  o 

(2.6)  B2  *  f  B2(y)dH(y)  +  f  [A(y)  -  A]2dH(y) 

Jo  Jo 

As  the  length  of  time  between  time  0  (when  the  event  TTq 

occurs)  and  the  time  that  the  next  particle  comes  and  produces 

an  impulse  is  distributed  as  F(x)  we  have 


(2.7) 

0(7)  -  M[ti  |  ^  =  ,3  -  H 

(2.8) 

D2(y)  »  D2p^|)(jj  *  j]  i.-2 

Theorem:  The  conditional  expectation  A(y)  «  |X^  *  y3  can 

be  determined  with  the  aid  of  the  following  integral  equation: 

(2.9)  A(y)  *f7  A(y-x)H(y-x)dF(x)  +f7  [  P  A(z)dH(z)3dF(x)  +  C(y) 
J0  J  0  y-x 


Proof: 

We  have 


"Pir*i  -% 


*  J,\  *  zj 


if  0  <  z  <  y-x  ahd  0<x<  j 
if  y-x  <  z  <oo  and  0«x<y 
if  y  <  x  <  oo  o 


and  (2.9)  follows  by  the  theorem  of  total  expectation. 


Theorem:  The  conditional  variance  B2(y)  «  IXq  =  can  be 

determined  from  the  integral  equation; 


(2«10)  B2(y)  *J^B2(y-x)H(y-x)dF(x)  +J*  [f°  B2(z)dH(z)3dP(x) 


+  [7[x+A(y-x)]2H(y-x)dF(x)  +  (7[JP°  [x  +  A(z)32dH(z)]dF(x)  +  D2(y) 
“0  ^  0  y «x 


Proof:  We  have 


D2!!*!  1^  -  Tq  *  *  y  j^i=z3 


B^(y-x)  if  0  <  z  <  y-x  and  0  <  x  <  y 
|B  (z)  if  y-x  <  z  <  oo  and  0  <  x  <  y 
.0  ify<x<oo. 


and  the  result  again  follows  by  the  theorem  of  total  expectation. 


It  should  be  recognized  that  equations  (2,9)  and  (2..10)  are 
of  the  same  type,  since  once  we've  solved  (2,9)  we  can  let 

P(y)  [x  +  A(y«x)32H(y~x)dF(x)  [x+A(z)J2dH(z)j  dP(x)  +  D2(y) 


for  this  will  simply  t©  a  function  of  y  „ 
Then 


o  „y  o  t>v  -  oci°  o 

BMy)  =1  Bc‘(y=>x)H(y-x)dF(x)  +j  {  j  Bw(z)dH(z)JdP(x)  +  r(y; 
J  0  J  0  y«x 


which  is  identical  in  form  with  (2,9), 


f 3«  The  application  of  the  method  to  special  dead  time  distributions. 

In  this  section  we  shall  treat  the  case  of  exponentially- 
distributed  dead  time  in  full  detail.  Then  we  proceed  to  formulate 
the  problem  for  max-m  impulse  times  and  solve  the  resulting  equations 
if  m=2  „  Finally  wo  show  how  the  Er'lang-2  case  leads  to  a  still 
mere  complicated  equation. 

Example:  H(x)  =  l»e“^iX 

We  let  Hfi-  *  e~ayA(y)dy  r(s)  »  J°  e“sydC (y) 

rr> 

c  f  e“aydF(y)  „ 

Jo 


and  fl(s) 


Equation  ( 2,9 )  becomes 


nJT  «<» 

(3ol)  A(y)  *f  A(y-x)H(y-x)dF(x)  +  C(y)  +  ’  I  A(z)dH(z)dP(x) 
JO  J0  y-x 


Taking  Laplace  transforms  wo  obtain 


s 


'i f  .  vir*  > 

-  V  /  v,  v«J>. 

sT  ^  ” 


A(s)  ^  P(gi,  ^  r-n/ 

S  3  *  S  ~  ^  ' 


i*  ' 


or 


(3.2)  'f  (s)  *-A(s)V(s+ij.)  +  A(s)M/{s)  +  H/UO-fL(s)  +  ru.) 


ar.  3  finally 


(3o3)  V  ;3*  A  =  ¥(A 


r>(s) 

‘ITTiT 


1*  *  A(s )  /  \ 

“  “7HFT  S  <!i) 


By  successive  applications  of  this  formula  v;e  can  express  Y's-i-rip.}  in 
terms  cf  YfA  and  ^(p.)0  This  lsad3  to 

(3oU  vr(:i+l'ii-l)ii)  “  ‘Piw  +  i  ~  i"il3+nu!  i  Vts+nii) 


if  v]  k.6  an  -  <tm  -  Jjfcgwfc  »r «  a=$$gj. 


this  simplifies  to  the  recursion  relation: 


(3o?)  +  (n+l)p)  =  d^  -  u>n  (|^(s+np) 


m  j  ic 


Ihe  /ven oral  solution  to  ( 3 0 5 )  j.< 


(3.,il  =  ia  <•  <-Uml'n  Vl  ... 


’rco.-r :  If  n  =  3.  tho  formula  (3,6)  satisfloa  (3,5)  .  /Isaume  Ws-Hrp) 

.  .ori  ::y  (3o  6J  for  k  « We  nov:  :»how  ^{a+Cn^Dp.)  is 


[V\  I,-  i;v 


f.?  '0.  %  the-  5  hrrt the tin 


-  <Wl  -h  <-l>  Vl  •••  %  ^=)  ♦  g 


0  o).  <i.T  (  —:l  J’1  J 


dubi  tltuliJm;:  ir  (3^0)  uo  obtain 


,  .V,.  (  ■ 

X 


cL^~  ■>  (“>!) 


wi  ••• 


a  (*q  Hi  3 


n-1  , , 

« ija+* 

•3RT  11  J- 


as  i  j  •:U‘oc",v 

V  f  '■  * 

,l  :i  *'•  t  a  J  i:i  a  Laplace  t  pansfoira 


liiu  *K(s+ntO 

rfy'OO 


-  0  «  Similarly  15m  jfKs+np)  =-  ■■) 


n-v'co 


Thus  lotting  a-> <m  In  (3*6)  i;o  conclude-  that 


Best 


Available  Copy 


If  we  set  a  =  n  we  find  that 


Thus  A{y)  can  be  determined  by  inversion  of  (3o7)o  It  is  of 
interest,  to  note  that 


(3o9) 


A  “  lo  A(y)y,e"tiydy  =  . 


The  conditional  variance  3  (y)  may  likewise  be  determined  by 
Laplace-Stielt jes  transforms 0 

Example  2:  H(x)  =  >  («1)^(?)  e"1^* 

J=0  J 

The  equation  to  be  solved  is 


py 

(3,10)  A(y)  =\  A(y-x)H(y«x)dF(x)+C  (y) 
J0 


A(z)dH(z)dF(x) 


Forming  Laplace-Stieltjes  transforms  we  obtain 


(3.1D  vi.)  =  g§}  +  ±(-dj9  [  twi- 

If  m  *  2  the  above  equation  becomes 


Y(3+ity.)  ss  2^(3+^)  +'H/'(s)<«0  +  Y(2p.)  -  2Y(yO  +  Yq 


where  „>  _  1— il( 

k - Ats-fgu) 


rk  s 


-ClsiMl 

TOiWT 


Theorem: 

"  [n”' J  fsum  of  products  of  exactly  k  ") 

(3.12)  YU+nn)  =  20-1-21'  fK®w1,...,»n_2  chosonj 

vouch  that  no  two  adjacent 
co^'s  are  taken 


Y(s) 


j~  n  -|  faum-s  of  products  of  exactly  k  ? 

~  2  n-2k  )terms  from  o  •  »^n-2  clloaen  \ 
'such  that  c*ij  ia  always  present/ 


/k- tuples  from  to^,.  ,d>n-2  such 

+  C  M^(2p)  *»  2V(p)3  (,  (2n*’^‘“l)  +  S_  \that  the  coefficient  of  the  l:-fcuple\ 

/ beginning  with  m  r  >  2  is  J 


•&1  °n-i<3+1) 


where  C  (s)  ** 
n 


{2n-2k-l[i.(  |  J1**1] 


^all  k- tuples  from  ,0 . .  ;«n<i>2  such  } 
gn-2k  Jthat  no  two  adjacent  4)^*3  are  chosen/ 
'form  products  and  add*  ' 


This  result  is  again  proved  by  induction.  The  inductive  argument 
is  omitted  as  this  approach  is  evidently  hopeless  and  we  did  not 
pursue  it  further. 

The  reader  may  find  it  interesting  that  the  corresponding 
equations  for  Erlang-r  dead  time  distributions  are  still  harder  to 
handle  than  the  ones  for  max-m  impulse  time. 

Example  3;  H(x)  =  1  -  ‘>™  e"^ 

the  basic  equation  (2,9) 

c7  P7  fifo 

A(y)  =  A(y-x)H(y-x)dF(x)  +  C(y)  +  \  |  A(z)dH(z)dF(x) 

J0  J0  fix 


now  becomes 


Taking  Laplace  transforms ,  recalling  that 

fl( s)  *  e"sxdP(x)  T(s)  =  F  e“aydA(y) 

OD 

and  that  n(s)  88  f  e“sydC(y) 

Jo 


we  obtain 
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£  1;.  Some  Final  Comments. 

From  the  practical  point  of  view  we  have  not  succeeded  In 
obtaining  useful  formulas,  however,  we  have  given  a  situation  where 
Erlang  distributions  are  definitely  not  to  be  used.  A  further 
interesting  observation  is  that  the  max-m  distribution  is  the  sum 
of  m  exponential  distributions  with  parameters  (m-l)p. ,. . . «  p.  „ 

Thus  a  distribution  which  is  the  sum  of  non-Identically  distributed 
random  variables  wa3  easier  to  handle  than  the  Erlang-distribution, 
the  sum  of  identically  distributed  random  variables  even  In  the 
Integral  equation  approach.  The  reason  for  the  success  of  the 
Erlang-distribution  in  the  theory  of  the  single  server  queue  is  that 
we  can  keep  track  of  the  "phase”  of  service  (or  arrival)  because 
the  customers  are  served  In  order c  Thus  the  the  behavior  of  E^/G/l 
can  be  deduced  from  that  of  the  single  server  queue  with  Poisson 
input  and  general  service  time  ’whore  the  service  Is  In  batches  of 
size  k.  In  the  infinitely  many  server  queue  the  order  of  the  arrivals 
Is  of  no  help  to  us  as  service  begins  on  arrival  and  the  second  arrival 
may  well  finish  his  service  while  the  first  is  still  in  the  queue. 

From  the  stand  point  of  trying  to  reduce  the  general  problom  to  one 
of  bulk  arrivals  the  max~m  distributions  would  seem  to  be  the 
distributions  to  use.  As  we  have  seen,  unfortunately ,  they  do  not 
simplify  the  problem  enough.  We  hope  to  have  some  more  to  say  about 
these  problems  at  some  future  data. 


Chapter  6 


A  Loss  Function  for  an  elementary 


missile  defense  system. 

In  a  talk  presented  at  Columbia  in  May  1961.  Professor  i?\rsek 

discussed  the  problem  of  defending  an  island  against  an  air  attack 
with  a  one  missile  battery 0  He  noted  that  the  system  has  an 
'absorbing  state”  in  the  sense  that  if  on  attacking  plane  is  not 

Or::yed  by  the  missile  battery  after*  tine  .  lo  will  be  divco'Xy 

over  the  base  and  will  attempt  to  destroy  the  base. 

In  thi3  note  we  shall  discuss  the  problem  of  determining  k'..r 
length  of  tine  until  absorption  (the  base  is  destroyed)  oe.-.  -.r:- 
assuming  that: 

1)  the  planes  arrive  according  to  a  recurrent  process,  I.e.,  the 

inter-arrival  times  tn+1  -  T  are  independent,  identically 
distributed  positive  random  variables*  ?far  ,  «.  'r  <  *'«  =  ■’!--)  i 

2)  The  time  it  takes  the  missile  battery  to  3hoot  a  plane  do on  Jo 

u  positive  raudou  variable  H(x)  „  If  "X  denotes  the  time  «t  takes 

II 

GO 

to  shoot  down  the  n«th  plane  then  ?  [Xn  «  xj  -  II U)  and  f$n  L-.}  are 

mutually  independent  and  independent  of  the  {Tn\  sequence,  be 
differ  from  the  usual  queueing  model  in  that  we  now  assume  that 
T[a>  t^J  *  1  ~  H(t0)  -  P  >  0  where  is  the  time  we  have  to  shoot 

the  plane  before  it  is  in  position  to  destroy  the  base* 

3)  We  assume  that  we  have  infinitely  many  guns  at  the  base  but  we 
use  only  one  gun  on  each  plans.  However,  we  always  have  a  giro 


This  approach  is  Just  that  of  taking  a  random  number  of  random 
variable 3 o  The  number  N  is  given  by  a  geometric  distribution. 

This  N  denotes  how  many  random  variables  ®n  **  ^  61,6 

chosen,  (see  Feller  £  X) )  • 

$2<>  A  slightly  generalized  Model 

In  the  previous  section  we  assumed  that  once  the  service  time 
(total  tine  used  to  shoot  the  plane)  reached  t q  the  plane  destroyed 
tho  base c  This  corresoonds  in  reality  to  perfect  accuracy  of  the 
bombardier©  Therefore  we  now  assume  bhat  when  tho  plane  is  over  the 
base  a  bomb  is  drooped  and  with  probability  C  >  0  It  destroys 
the  base<»  In  order  to  determine  Cr(x)  in  thi3  case  we  simply 
redefine  p  equal  [l  -  H(t0V}c  which  now  gives  the  probability 

that  the  n-th  plane  destroy r-  the  system. 


Theorem  2; 


00 

Jo 


-8X.,. 


du(x) 


C(l  -  H(t0)]  e"sfc0 

i  -  c€(s)Ci-H(t0)] 


A  decision  theoretic  criteria  for  efficiency. 

In  statistics  one  criteria  for  deciding  upon  a  good  estimate 

2 

is  to  choose  a  los3  function  (usually  squared  error  L(6,@)  »  (5(x)-8)U 

and  using  that  estimate  minimizing  the  expected  loss.  Unfortunately 
in  most  queueing  problems  it  is  very  difficult  to  calculate  these 
expected  losses.  However,  for  our  model  it  is  easy.  ^  fairly  realistic 
loss  function  would  be  of  the  form  T.  s  Loss  =  Wf  *  p(t'  where  k 
would  be  a  constant  such  as  cost  per  plane  and  p(t)  would  be  a 
polynomial  in  t(tima). 


E(L)  *  kB(N)  +  E(p(t) 


*  ■  0 

But  then 


75 


*  ~  +  E(p(t) ) 


But  from  the  fact  that 


Jo 


e“8XdO(x)  = 


we  can  calculate  all  the  moments  of  the  random  variable  t  thus 

can  calculate  E(p(t)),  For  instance  if  p(t)  *  t 
t~p  +  q*< 

E(t )  is  -£-* - and  thus  if  L  =  kN  +  t 

P 


E(L) 


kp  +  t^p+  q*< 


It  is  easily  seen  that  the  k-th  moments  of  G(x)  depend  only  on 
p,  tQ  and  the  first  k  moments  of  F{x)  .  Thus  If  all  moments  of 
F(x)  exist  so  do  all  the  moments  of  0(x)o 
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